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ABSTRACT 

Giant star-formation clumps in dwarf irregular galaxies can have masses ex- 
ceeding a few percent of the galaxy mass enclosed inside their orbital radii. They 
can produce sufficient torques on dark matter halo particles, halo stars, and the 
surrounding disk to lose their angular momentum and spiral into the central re- 
gion in 1 Gyr. Pairs of giant clumps with similarly large relative masses can 
interact and exchange angular momentum to the same degree. The result of this 
angular momentum loss is a growing central concentration of old stars, gas, and 
star formation that can produce a long-lived starburst in the inner region, identi- 
fied with the BCD phase. This central concentration is proposed to be analogous 
to the bulge in a young spiral galaxy. Observations of star complexes in five 
local BCDs confirm the relatively large clump masses that are expected for this 
process. The observed clumps also seem to contain old field stars, even after 
background light subtraction, in which case the clumps may be long-lived. The 
two examples with clumps closest to the center have the largest relative clump 
masses and the greatest contributions from old stars. An additional indication 
that the dense central regions of BCDs are like bulges is the high ratio of the 
inner disk scale height to the scale length, which is comparable to 1 for four of 
the galaxies. 
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Subject headings: Galaxies: bulges — Galaxies: evolution — galaxies: irregular 
- Galaxies: starburst — galaxies: individual (Mrk 178, DDO 155, Haro 29, 
NGC 2366, NGC 4861) 



1. Introduction 



Blue Compact Dwarfs (BCDs) are sm all galaxies with intense emission- lines from star- 
burst HII regions in t heir central regions (Sargent fc Searle Il970 ). They are gas-rich like 



other dwarf irregulars ( Chamaraux 1977 : Gordon fc Gottesmanl 



but m uch more centr ally concentrated in stars (Noeske et a. 



2006), star formation (iHeller et al.ll2000t iHunter &: Elmegreen 



981 



Thuan fc Martin 



2003; 



1981) 



Hunter fc Elmegreen 



20041 ). and gas (ITavlor et al. 



1994; van Zee et al.lll998bh . This concentration suggests that gas inflow following angu- 



lar momentu m loss led to enhanced star formation in a dense and gravitat ionally unstable 
central disk ( ITavlor et al.lll994l ; iPapaderos et al.lll996t Ivan Zee et al.l 120011 ) . Gas loss from 
the outer parts also produces a shrink ing radius for star formation, as observed in dwarfs 
rtzhang et allhoill iKoleva et al.lboilh . 



Individual star-forming regions in BCDs are relatively large, giving the galaxies a 



clumpy, irregular a ppearance in Ha (e.g. iKunth et al.l Il988t ICairos et al.l l2009bl ) and FUV 



1981; 



Firpo et al. 



( Th uan et al.lll997 ). Clu mp emission lines are supersonic and apparently virialized ( ITerlevich fc Melnick 



201 ll ) . w hich implies the clumps could last for several internal cro ssing 



times. Low velocity shear dThuan et al 



1999; 



van Zee et al.l l200ll : iRamya et al.l l201ll ) and 



resolved stellar population studie s (Dohm- Palmer et al. 19981) also suggest the c lumps could 



be long-lived, 100 Myr or more. McQuinn et al ( 201ol ) and Zhang et al. ( 2011 ) suggest the 
starburst itself can last for ~ 1 Gyr. 

S urrounding many BCDs are pools of HI , sometimes as large as 4 or more optical 



radii dBrinks fc Klein 



Pustilnik et al. 



ing (Ivan Zee et al 



1988; 



Taylor et al. 



1996; van Zee et al. 1998b: Putman et al. 



1998 



20011: iHoffman et al. 

1998bh or even in a disk 



20031 ). The peripheral 



Pustilnik et al 



TI is often not simply rotat- 
1997I. boOlh . This peripheral 



gas led to suggestions about cloud impacts dGordon fc Gottesmanl 1 1 98 ll) , wea k interaction s 



dBrinks fc Klefn]|l988[ IPustilnik et aDbopj iBravo-Alfaro et al.ll2004l) . merging jBekkill2008h 



and tidal effects (Ivan Zee et al.l Il998bt IPustilnik et al.ll200ll ) in efforts to explain the high 
star forma tion rates. Alterna t ively, the HI could be vestigial streams or po o ls of cosmological 



accre tion ( ITavlor et al.l Il993t iThuan fc Izotovl Il997t Ivan Zee et al.lll998at IWi 



19981 ). like the streams modeled for higher-mass galaxies (e.g., 



Ceverino et al. 



cots fc Miller 



2010|). 
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BCDs resemble youn g; galaxies in many respects (lElmegreen et al.ll2 009b 



2011: 


Griffith et al. 


2011) 


turbulent ( 


van Zee et al. 



1998b. 2001 



200 



Silich et al. 2002; Garcia-Lorenzo et al. 



zotov et al. 



), rela tively 
2008h . and 



highly clumped with locally intense star formation. They tend to lie at the edges of galac- 
tic clusters or in voids with only low-mass neighbors, and where harsh environmental ef- 



fects like ram pressur e stripping are minimal ( iGrogin fc Gellerll2000l ; iDrinkwater et al.ll2001 



Pustilnik et al.ll2002l ). Those with the lowest metallici ty tend to be relatively young, having 
formed most of their sta r s within the last few Gyrs (ISearle fc Sargentl 1 19721 ; iNoeske et al. 



20001 ; IJohnson et al.ll2000l ; iFricke et al.ll200ll : IPapaderos et al 



metallicities tend to have relatively more old stars (e . 



Crone et al. 2002; Caon et al.l2005: Cairos et al.l2007 



Thuan 



2009a 



2008). BCDs with less extreme 



1983; 



Loose fc Thuan 



1986 



b uZhao et al.ll2011uZhang et al. 
20111 ). In the most low-metal BCD s, there is little evidence for stars older t han ~ 200 — 500 



Myr; these include S BS 1415+437 ( 



jGuseva et alJl2003bh . and I Zw 18 JPapaderos et al.ll2Q02h 



huan et al.lll999l: iGuseva et alJl2003ah . SBS 1129+576 



Specific examples of BCDs illustrate these points. I Zw 18 is in many respects a mor- 
phologically young system. It has two giant star-forming regions inside a kpc-scale blue 
conti nuum of stars, ionized emission with an overall exp onential profile (IPapaderos et al. 
20021 ) . and an extensive HI e nvelope (Ivan Zee et al.l Il998al ) . Inside each region the star for- 
mation is widely distributed (IHunt et al.l 120051 ) . The rotation curve is flat but steeply rising 



201 lh. Ra dial 



mo- 



in the inner part, where baryons may dominate dark matter (ILelli et al 
tions in the disk of ~ 15 km s _1 suggest a major disturbance, which Lelli et al. ((2011) suggest 
is tidal because there is a dwarf companion galaxy and elongated peripheral HI gas. The 
rotat ion curve gives a mass of 10 8 M® in which ~ 70% is neutral gas (IContreras Ramos et al. 
201ll ). It h as very low metallic ity (2%-3% solar, llzotov et al.l 120011 ) even though old stars 
are present. iRecchi et al.l (120041 ) and others have considered the selective removal of metals 
by winds. 

Another well-studied example is VII Zw 403 , whic h is among the nearest BCDs to the 
sun, having a distance of 4.5 Mpc (ILynds et al.l 119981). VII Zw 403 has a half-dozen big 



a smooth elliptical background of old red giant stars dSchulte-Ladbeck et al 



clumps and many Ha filaments from superbubbles (lLozinskava et al.l 120061) in the m idst of 

The 
The 



dynamical mass is 2 x 10 s M , with a pproximately 20% in HI (IThuan &; Martin 



998). 



19811 ) 



rotation speed is only ~ 15 km s " 1 (ISimpson et al.ll201ll) and the metallicity is 5% solar 
(Martin 1997 ; Izotov et al. 1997 ). Lynds et al. (jl998 ) used resolved stellar populations to 
date a major star burst to 600 Myr ago, when the star formation rate was ~ 30 times higher 
than it is today. The current burst produced ~ 10 6 M within the last 10 Myr (jSilich et al. 
2002h . 
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The most intense star formation in BCDs can occupy very compact regions with extreme 
densities and lo cal formation rates. SBS 0335-052 is a pair of extremely young interacting 
dwarf galax ies (Pustilnik et al.ll200ll ; lEkta et al.l 120091 ) without much of an u nderlying old 



population (IPapaderos et al 



2009a; Peimbert fc Peimbert 



1998h . and with a metallicity of 2.5%- 4% solar (llzotov et al.l 



201dh . There are 6 super star clusters ( jThuan et al.lll997h. of 



which two, within ~ 200 p c of ea ch other, have extremely intense star formation. iHunt et al. 



( 120051 ) and iJohnson et al.l ( 120091 ) found radio free-free absorption and a very high emission 
measure where the electron density is ~ 10 3 -10 4 cm -2 , the star formation rate is ~ 1 M F 
yr _1 , and the excitation comes from the equivalent of ~ 10 4 07 stars. 



Here we propose that central accretion and long-lived starbursts in some BCDs arise 
from gravity-driven motions and torques produced by clump formation, clump dynamical 
fricti on, and clump interactions - the same processes that coul d make bulges in larger galaxies 
(e.g. lNoguchilll999l ; llmmeli et al.ll2004l ; iBournaud et al.ll2007f). BCDs hav e steep stellar pro- 
files in the inner 500 pc that are exponential (IHunter fc Elmegreenll2006l ) or deVaucouleur's 
(IDoublier et al.lll999l ). as in the bulges of earlier Hubble types. Such high central concentra- 
tions require baryonic mass inflow and significant angular momentum redistribution in the 
disk. Much of this inflow could have occurred when BCDs were young, but some of today's 
BCDs still look dynamically young even if there are old stars, and significant inflow could 
be occurring now. 



Our emphasis differs from that in iGovernato et al.l (120101 ). where simulations of dwarf 
galaxies highlight the removal of gas in order to avoid central concentrations. In these sim- 



ulati ons, diffusion, tore 



ues, and pressure-driven inflow s return some of this gas to the center 



Recchi fc Henslerl 120061 ; iDalcanton fc Stilpl 12010 ) , only to haye it removed again by th e 



(e.g. 

next starburst, cycling in and out many times (e.g. 
The result is a bulge-free late-type galaxy and a time-changing centra 
a primordial dark matter cusp into a more uniform dark matt er core ( 



Stinson etal.1 120071 : iRevaz et all 120091 ) . 



The de gree of t his co nversion varies for different simulations (jOgiya &: Mori 



potential that convert s 
Read fe GilmoT]l2005h 



2011 



Oh. et al. 



201ll ). lAlardl (120111 ) note that the least evolved galaxies, having the highest relative gas 
abundances, tend to have the steepest inner density profiles, supporting the idea that gas 
recycling and star formation make the inner profiles shallow over time. These observations 
could imply that so me BCDs sti l l have steep central dark matter profiles, if these galaxies 
are relatively young. iDel Popolol (120111 ) also model low-mass galaxies and suggest that tidal 
torques and the baryon fraction before collapse influence the central density profile of dark 
matter. Steeper central profiles are predicted to occur in more remote gala xies and in those 
with higher dark matter fractions; BCDs could be in this category too (IGrogin fc Geller 
2000h . 
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BCDs are unusual in having both a central concentration and a hig h gas abundance. 
This c ombination also appeared in massive galaxies at redshift z ~ 2 (jElmegreen et al. 
2009aj). Observatio ns at intermediate-tq - high redshift indic ate that Hubble types arise 



mostly since z ~ 1 (IPapovich et al.l l2005r iBundy et al.ll2006l ). The morphology seems to 



depend on dynamical maturity. High- mass dis ks like Hubble type Sa ten d to be higher den- 
sity than low- mass disks like Hubble type Sd ( jRoberts fc HaynesHl994j ). Thus, Sa's evolve 
more quickly to a c entrally concentrated state with a low gas fraction and a low specific 
star formation rate ( jSandagd Il986l ; IZhang fc Butal 120071 ) . If this trend of increasing central 
concentration and decreasing specific star formation rate continues into the future for low- 
mass galaxies, then some dwarf Irregulars might also evolve to a centrally concentrated state 
with little remaining gas. Mass loss from winds and supernovae in the low potential well 
of the dwarf could prevent such a central concentration however, depending on the relative 
rates of inflow from torques and outflow from winds. Those with dominant inflows could 
go through the BCD phase, as discussed here. Strong inflow depends on the presence of 
relatively massive clumps or tidal arms. 

In what follows, we estimate the accretion time of a clumpy disk from dynamical friction 
(Sect. 2), and then consider whether observed clumps and other irregularities in BCDs are 
massive enough to drive significant disk evolution on a Gyr time scale (Sect. 3). A summary 
is in Section 4. 



2. Clump Accretion in BCDs 
2.1. The Case with a few Giant Clumps 

The process of clump drag and interaction leading to coalescence i n the center of a galaxy 



has been illustrated w ith detailed simulations of high redshift galaxies (jElmegreen et al.ll2008 



Ceverino et al.ll2010l ). The clumps in these simulations formed spontaneously in a turbulent 



disk and had masses of about 5% of the total galaxy mass. Migration to the cent er took only 



a few orbit times (~ 0.5 Gyr). Clump destruction by star formation feedback (IGenel et al. 



20101 ) does not stop the torques and accretion if each destroyed clump is replaced by a new 
one. This replacement is likely as long as the conditions for forming the first clumps, such 
as high gas fractions and turbulent speeds, are still present. A low ratio of turbulent speed 
to orbit speed would stop this process, because then the clumps that form by gravitational 
instabilities are relatively low-mass and produce proportionally weaker torques. The biggest 
star-forming regions in the Milky Way are only ~ 10 -3 times the disk mass and should have 
little tendency to move to the center. 



-6 - 



Dynamical friction and clump torques are important if the ratio of the disk Jeans mass 
to the galaxy mass is more than a few percent. This ratio scales with the square of the ratio 
of the gas velocity dispersion to the rotation speed. In high redshift galax ies, the rotation 



speed is normal for a massiv e disk but the dispersion is abnormally high (jErb et al.l 12006 



Forster Schreiber et al.l 1201 ll . and references therein) making the ratio high. In local dwarf 
Irregulars, the gas dispersion is normal for local galaxies, ~ 10 km s _1 , but the rotation speed 
is low, ~ 50 km s _1 or less. In both cases the ratio of speeds is high and the clumps that form 
by gravitational instabilities are massive compared to the disk. The same processes of massive 
clump formation and angular momentum exchange should happen in high redshift galaxies 
and local gassy dwarfs because both have relatively large velocity dispersions compared to 
rotation speeds. 

The timescale for dynamical friction between an orbiting clump of mass M c in the disk 
and non-rotating dark matter particles or stars in a halo is v (dv/dt) -1 where 



dv 
~di 



47rlnAG 2 M c p 



erf [X] 



2X 

—2< 



-X 



71 



4vrlnAG' 2 M c p£ 



( Binney fc Tremainen2008l ). Here, v is the clump orbital speed, X = v/(2 1 ^ 2 a) for halo 3D 
velocity dispersion a, p is the halo density, In A is the coulomb factor, and £ is the quantity 
in parentheses. This formula assumes that the clump is a self-gravitating object surrounded 
by a uniform density of low-mass field stars or dark matter particles that have a Maxwellian 
velocity distribution function. 

It is convenient from an observational point of view to write the local rotation speed as 
a power of the local radius, v(r) oc r 13 , since /3 comes from the rotation curve. Starting with 
p(r) = p or - a and v(r) 2 = GM dyn (r)/r, we get v(rf /(Anp[r]Gr 2 ) = 1/(3 - a) = 1/(1 + 2(3). 
Then the dynamical friction time, v/(dv/dt), in units of the dynamical orbit time, r/v, is 



T(r) 



M, 



dynl 



r(dv/dt) In Af (1 + 2/3) 



M r 



T n (r) 



M, 



dyn\ 



M r 



(2) 



where Md yn is the galaxy dynamical mass enclosed within the orbital radius of the clump. 



Dwarf galaxies have nearly solid body rotation in the inner parts (ISwaters et al.ll2002l ). 
BCDs can have steeply ri sing rotation curves in the dense inner regions, and flatter rotation 
curves beyond that (e.g., Ivan Zee et al.lll998bl ; iLelli et al.l 1201 ll ). The BCDs we consider in 
Section 3 have approximately-linear rising rotation curves in the vi cinity of the giant c lumps, 
and some have flat rotation curves beyond that (e.g., NGC 2366; iThuan et al.ll2004l ). Thus 
for the main starburst regions we can take /3 ~ 1 or slightly less. If v ~ a, then £ ~ 0.20. In 
that case, the dimensionless time coefficient in equation ([2]) is To = 0.56 for typical In A = 3 
(see below). For a clump at r = 0.5 kpc orbiting with v — 10 km s" 1 , r/v — 49 Myr and the 
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dynamical friction time in physical units is 27 'M^ ya /M c Myr. This means that a clump with 
a mass greater than 2.7 percent of the enclosed galaxy mass has T < 1 Gyr. 

Note that Md yn decreases with radius, while the normalization quantity, r/v, is constant 
for /3 ~ 1. Thus the timescale gets smaller as the clump moves in. Writing the rate of 
change of clump angular momentum as dL/dt = M c r(dv/dt) for frictional deceleration in 
the azimuthal direction dv/dt, and setting this equal to M c v(dr/dt) for circular speed v and 
radial drift speed dr/dt << v, we get dr/dt = v/T(r). For /3 = 1, v oc r and Md yn oc r 3 so 
T(r) oc r 3 if A and £ are constant. Then, the time to reach the center is 1/3 the instantaneous 
T in equation [2j If f3 — 1/2, then v oc r 1//2 and Md yn oc r 2 so dr/dt oc r _3//2 and it takes 0.4T 
to reach the center. 

One uncertainty in this result is the ratio of the disk orbit speed to the halo velocity 
dispersion, which enters into £. This ratio depends on whether t he halo has a c ore or a cusp, 



and on the nature of the core. We consider t wo extreme cases: alBurkertl (119951 ) halo density 



profile in the case of a core, and an NFW (jNavarro et al.lll996l ) profile for a cusp. For the 



Burkert profile, p(x) = p s ([1 + x][l + x 2 }) 1 with x = r/r s and scale factors p s and r s . We 
use this with the equation of hydrostatic equilibrium in the radial direction to determine 
(3 and v/cr as functions of x. Hydrostatic equilibrium implies dP/dr = —GM(r)p(r)/r 2 
where P = pcr\ D for ID dispersion o\d = a/3 1 / 2 and M(r) = J Q r Amr 2 p{r~)dr . We assume the 
boundary condition P ~ and a ^constant at the edge of the halo, which is taken to be 
where p = 10~ 4 po- Figure [TJ shows v/cr, £, (3, and T as functions of position x (determined 
by numerical integration). ~ 1 for a solid body rotation, so x must be small in the visible 
part of the disk. For example, x = 0.54 at the half-density point, where p = 0.5po, and there 
v/a = 0.49, £ ~ 0.03, (3 = 0.72, and T ~ 4.6. For r = 0.5 kpc and v ~ 10 km s" 1 , the 
friction time is then rT /v = 220Md yn /M c Myr. This implies that a clump with a relative 
mass of M c /Md yn = 5% takes ~ 1.5 Gyr to spiral in, considering the factor of 1/3 that 
accounts for a decreasing M dyn with radius, as discussed above. 

The NFW dark matter profile is a little faster. For this, p = p s (x[l + x] 2 ) 1 with x = 
r/r s again. There is a density singularity at the center that produces a logarithmic divergence 
of a quantity like pressure if the equation of hydrostatic equilibrium is considered (because 
M(r)pjr 2 ~ 1/r near the center). We consider instead tha t the halo 3D velocity dispersion 



is comparable to the rotation speed (jNavarro et al.l Il9961 ). v ~ a, which gives £ = 0.20. 



Writing for the galaxy mass M^ yQ (x) = 4nr^p s A4(x), where A4(x) = ln(l + x) — x/(l + x), 
the slope of the rotation curve is now given by 1 + 2f3 = x 2 / (A4(x)(l + x) 2 ). Figure [1] shows 
(3 as a dashed line. We see f3 ~ 0.5 for small x and then T ~ 0.83 with £ = 0.2. This is a 
smaller prefactor than for a cored halo because of the higher v/cr. The dynamical friction 
time for an NFW halo is ~ 41Md yn /M c Myr with r = 0.5 kpc and v ~ 10 km s _1 . Clumps 
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with M c /Md yn > 1.6% take < 1 Gyr to spiral in, considering the factor of 2/5 that accounts 
for a decreasing M dyn with radius when (3 = 0.5. 



Dwarf galaxies do not appear to have NFW halos at the present time (e.g., 


Zackrisson et al. 


2006; 


van Evmeren et al. 


2009; 


Kuzio de 


^arav 


2009; 


Amorisco & Evans 


20121), although it 


is difficult to be certain 


(Spekkens et al. 


2005: 


Valenzuela 


2007 


) and it is not yet known 



whether BCDs differ from other dwarfs in this regard. An important point for dynamical 
friction is the relative velocity dispersion of the field particles, cr/V, which en ters into £ as 



shown above. The validity of the Chandresekhar formula is also a question. iGoerdt et al. 



( I2010l ) modeled sinking massive objects like what we consider here and showed that a cen- 
tral cusp turns into a core inside the radius where the sinking mass equals the enclosed dark 
halo mass. Further sinking in their model stalled at this rad ius because of a d ecrease in 
dynamical friction in the core. This result was also found by iRead et al.l (120061) and oth- 
ers. In an detailed study of dynamical friction in cored galaxies, Inouel ( 2011 ) explained the 
loss of frictional forces as a result of orbit resonances that appear when the orbit time is 
independent of radius, as is the case for a constant central density. Another limitation is 
that after a bulge forms, tidal forces fr om the bulge can rip apart remaining clumps and 
prevent them from reaching the center (jElmegreen et al.l 120081 ). These considerations make 
it plausible that in some BCDs today, giant clumps come in from larger radii and then stall 
at the edge of a dark matter core, building up the central region from inside out. Accretion 
to the center could have occurred more readily when the galaxies were young and the dark 
matter profiles were more cuspy. 



Sanchez-Salcedo et al.l ( 120061 ) considered a similar situation with in-spiraling globular 
clusters in dwarf galaxies. They concluded that dwarfs should have dark matter cores rather 
than cusps so that the globular clusters stall midway in the disk where they can still be seen. 
Galaxies with gi ant nuclear s tar clusters, however, may have allowed th eir disk clusters to 
reach the center (lBokerll2010l ; lAntonini et al.l 1201 It lHartmann et al.ll201ll ). The resemblance 
of these galaxies to BCDs is compelling if we allow for a difference in the mass and size of 
the disk clumps that form by star formation: BCDs have relatively massive and large disk 
clumps that could spiral in to form massive and large central cores, while normal dwarfs 
and galaxies have relatively small star clusters, which could spiral in to form small nuclear 
clusters. 



2.2. Clumps with a Power-Law Mass Function 



As mentioned above, individual clumps need not survive the full trip to the center to 
drive accretion if new clumps replace dispersed clumps in a steady state. Similarly, there 
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need not even be a single giant clump. Any irregularities moving through a slower rotating 
halo will have dynamical friction drag, and the total torque on the medium will depend 
on the mean squared mass of those irregularities. Consi der the equat ions of disk accretion 
starting with the continuity equation in two dimensions (jPringle 1981 ): 



d -i d _ , , 

-S + r" 1 — rS^ r = 0; 3 
at dr 

the surface density is E and the radial drift velocity is v T . The torque equation is 

-^-rEt> + r -1 -^-rJ2v T rv = torque/area (4) 
dt dr 

where v is again the azimuthal speed. For a viscous disk, the torque per unit area is 
R~ x dG/dR where G = uHAr 2 for viscous coefficient v (comparable to the product of the 
clump mean free path and the rms speed), and Oort rotation constant A, which is the rate 
of shear. We are not concerned with viscosity in this paper because A is small for dwarf 
galaxies with little shear; v is also u sually small compared to dynamical torques in galaxy 



disks (however, see lWang et al.ll2009l ). Here we consider dynamical friction between orbiting 
clumps in the disk and the halo and other parts of the disk. If an annular area has iV 
clumps of mass M c , then the total torque exerted on these clumps is N M c rdv / dt for decel- 
eration in the azimuthal direction dv/dt from equation (CQ). This quantity depends on the 
mean squared clump mass because dv/dt depends on clump mass. For clump mass function 
dn(M c )/dM c cx M~ s , 

2 6 

< NM 2 >= —— M c , max < NM C >= /M c , max < NM C > . (5) 
6 — o 

The prefactor f range s between / = 0.33 at 5 = 1.5 and / = 0.06 at 5 = 2 (for 5, see e.g. 
Heithausen et al.lll998l ). In the latter case, the integral over M 2 n(M) gives ln(M Cjmax /M C)r 



:,mm 1 1 



whose value is ~ 16 for typical M c . 

max r ^ J 

10 7 M and M c . 

min n * J 

1 M Q . Because < NM C > 

/Area = E, we have < NM 2 > /Area = /EM c max . Then the torque/area for the above halo 
model becomes fErdv/dt = Ej/r where 7 = (1 + 2/3)/ In AGM cmax ^. For a fixed galactic 
potential {y independent of time), the torque equation is now 

r 2 v— — h — rEt> r rt> = — E7. (6) 
dt dr 

We can simplify this by writing fi = rEi> r and noting that (d/dr)firv = fi(d/dr)rv + 
rvdfi/dr = n(d/dr)rv — r 2 v(d'E/dt) using the continuity equation. Then dY^/dt cancels 
in the torque equation and we get 



(7) 
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This may be solved for v z since everything else is a known function of r, and then the result 
can be put into the continuity equation to get dYi/dt. Note that drv/dr = (1 + f3)rv; for the 
other radial derivative, we set d/dr ~ 1/r. The resultant normalized accretion timescale is 

T _ (1 + /3) M dyn (x) , M dyn (x) 

rc/S/rft ~ /lnA£(l + 2/3) M c , max U J M c , max ' 1 ] 

This time is larger than before by the ratio (1 + (3)/f. Setting = 1, In A = 3 and £ ~ 0.03 
for a Burkert core, and taking / ~ 0.1, we get T\ = 74 and a physical accretion time of 
1.2Md yn /M Cjmax Gyr for r = 0.5 kpc, v — 10 km s" 1 with the factor 1/3 to account for a 
decreasing Md yn with radius. For a NFW core with = 0.5 and £ = 0.2, T\ = 12.5 and 
the accretion time is 240Md yn /M c max Myr for r = 0.5 kpc, v — 10 km s _1 with the factor 
2/5. Now we see that it takes about 1 Gyr for 10% of the ISM to accrete to the center from 
the inner half-kpc if the largest cloud in a power law distribution of cloud masses is 12% 
and 2.4% of the enclosed galaxy mass for the Burkert and NFW profiles, respectively. If 
the largest cloud has a mass much larger than the extrapolation of a power law distribution 
from the other clouds, then the previous analysis for a single cloud applies. 



2.3. Clump-Clump Interactions 



Clump-clump interactions can also drive accretion by direct gravitational forces. The 
acceleration on one clump by another clump is GM C /Ar 2 for separation Ar. This acceleration 
cumulatively distorts the clump's motion until its velocity has changed significantly. The 
timescale for this change is v divided by the acceleration, and in units of the orbit time, it is 
v 2 /r divided by the acceleration. Writing v 2 /r = GM^ yn /r 2 , the normalized interaction time 
becomes (Ar/r) 2 (Md yn /M c ). For big clumps, Ar ~ r, and the normalized interaction time 
is just the ratio of masses. That means To ~ 1 in an equation like ()2]), and the accretion 
time is ~ 49Md yn /M c Myr for r = 0.5 kpc and v ~ 10 km s -1 . 

Other torques will contribute to the inflow of gas, so T is an upper limit based on 
dynamical friction with the halo. The disk als o will produce a torqu e if there is shear, 



because then the clump will drive a spiral wake (Juli an fc Toomre 



1966) and the wake will 



have its own torque that drives mass inward ( jLynden-Bell fc Kalnajslll972l ). Small galaxies 
tend to have little shear, however. 
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3. Observations of Clump Properties in Dwarf Irregulars and BCD Galaxies 
3.1. Clump Mass Fractions and Accretion Times 



Many of the observations referenced in the introduction concern BCD or other dwarf 
irregular galaxies that have relatively large gas velocity dispersions compared to the rotation 
speeds, relatively thick disks compared to the galactic radii, and relatively massive star- 
forming regions compared to the galaxy masses. They are good candidates for the extreme 
torques and inward migrations discussed above. Five examples are given in Table [H along 
with properties of their primary star-formation clumps. As above, the clump mass is denoted 
by M c , the galactocentric radius at the center of the clump is r c , and the galaxy mass 
inside the clump radius is Md yn , measured as r c v(r c ) 2 /G for local rotation speed v(r c ). Also 
for reference, we give the total galaxy stellar mass, M s , the total baryonic galaxy mass 
(gas+stars), M b , the Coulomb factor A ~ {R dump / R z ^ ot )(R g ^ to Jr c ) l+2 ^(M Ajn / M c ), and 
the accretion time T from equation (j2J) for NFW and Burkert dark matter profiles. The 
Coulomb factor comes from the approximate expression A ~ (-Rciump/-Rgai,tot)(-^gai,tot/-^c) in 
Binney &: Tremaind (120081 ). where -R c i ump and M c are the clump radius and mass, and i? ga i,tot 
and Mga^tot are the total galaxy radius and mass. We take M ga i itot /Md yn ~ (-Rgai.tot/^c) 1+2/3 
for rotation curve v oc r 13 , -R c i ump equal to half the clump aperture in the table, and -R ga i,tot/ r c 
from Figure [2j As discussed above, the NFW profile has j3 = 0.5, £ = 0.2, and a time 
multiplier of 0.4 to account for faster accretion as the clump approaches the center; the 
Burkert profile has (3 = 1, £ = 0.03, and a multiplier of 0.33. 

Color composite im ages are shown in Figure El F or four galaxies, they are made with U 
and J band images from lHunter &: Elmegreenl (120061 ); for NGC 4861, they are made with u 
and z band images from the Sloan Digital Sky Survey (SDSS; Stoughton et al. 2002). The 
elliptical contours outline the giant star-forming clumps that we consider to be candidates 
for inward migration. The cross marks the center of the galaxy as defined by the outer 
elliptical isophotes in V band. 

The HI rotation curves for our galaxies are close to solid body in the clump region. 
They are not accurate enough to tell if the dark matter halo is cored or cuspy (which give 
rotation curve slopes j5 — 1 or 0.5 in these two cases, respectively). To determine v at the 
clump radius r c , we fit the observed rotation curve to a deprojected speed, v t , and radius, 
r t , at the limit of the observation or the turnover point, whichever comes first. Then the 
total galaxy mass inside the clump radius is taken to be v 2 r^/(r 2 G), assuming = 1 for 
this. For the different galax ies, the values of (v t : r t ) in (km s" 1 , k pc), are, DDO 155: (8, 
0.151 JCarignan et al.l[l99oh: Haro 29 : (34, l^jstTfebrJ ^S NGC 2366 : (50, 3.3) 



(IHunter et al 



2001 



Thuan et al.ll2004 ): NGC 4861: (40, 3.3) (IThuan et al.ll2004h . 
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For DDO 155 (also known as GR8), the HI observations by lBegum fc Chengalurl (120031 ) 
suggest t h at th e velo city field is complex , so th e dynamical galaxy mass is inaccurate. 
Lo. et al.l (Il993h and iBegum &: Chengalur fl2003h suggest the y elocities have an expand- 
ing or contracting component, which IBegum &: Chengalurl (120031 ) fits to a peak value of 10 
km s -1 , with some radial variations. They also fit the rotating part to a peak value of 6 km 
s _1 at the edge (the escape speed was estimated to be ~ 30 km s^ 1 ). The velocity of the 
giant clump outlined in Figure [2] is smaller than the systematic velocity, so if it is on the near 
side of the galaxy, then it is expanding away from the center. IBegum fc Chengalurl (120031 ) 
consider an explosive origin for this motion but note the lack of old star clusters that might 
have driven this explosion; the y suggest that HII reg i ons m ight have had the necessary force. 
If the motion is inward, then IBegum fc Chengalurl ( 120031 ) suggest that the clumps might 
have coalesced to form the galaxy and are now dispersing to make a disk. They note that 
there are no tidal features, however. This inward moving interpretation is consistent with 
the model presented in the present paper; tidal features are not expected because the torques 
are generated internally. The timescale for inward motion giv en in Table 1 is ~ 2 My r 
for DDO 155. This timescale is consistent with the results of iDohm-Palmer et al.l ( 119981 ). 
who find from resolved stellar population studies that star formation lasts in each clump 
for ~ 100 Myr. They suggest that this long time requires gravitational self-binding of the 
clumps. Dohm-Palmer et al. also suggest that the clumps come and go on this time scale, 
with the current generation of clumps at the positions of the three main HI clouds. Thus the 
present model of massive clump formation by gravitational instabilities in a gas-rich galaxy, 
relatively long clump ages from modest gravitational self-binding, and angular momentum 
loss through halo, disk, and clump-interaction torques, is consistent with the HI and stellar 
observations and previous interpretations of DDO 155. 

Mrk 178 does not have a published rotation curve, but the HI line width was given by 
Bottinelli et al.l (119731 ). who also derived a total dynamical mass ("indicative mass") from 
the equati on M tnt civn = 3 x 10 4 ? -hH /2 M q for Holmberg radius r H in kpc and linewidth W 
in km s _1 (IBottinelli et al.lll968l ). Scaling to our distance, the Holmberg radius is rn = 2.9 
kpc and the total dynamical mass is Mh = 1.6 x 10 9 M & . If we assume this rotation curve 
is solid body, then the dynamical mass inside radius r c is M dyn (r c ) = M H (r c /r H ) 3 . Setting 
r c = 0.39 kpc from Table 1, we get M dyn (r c ) = 3.9 x 10 6 M . 

Clump and galaxy stellar masses were derive d by fitting the S EDs over a range of 
passbands inside deprojected circular apertures (see lZhang et al.l 1201 ll ) . The aperture sizes 
were determined from the U-band brightness contours shown in Figure 2. Local background 
intensities came from larger annuli around the clumps and were subtracted from the clump 
intensities. For M rk 178, DDO 155 and NGC 2366, the SEDs used observations in U, B, V, 
and J passbands (IHunter fc Elmegreenll2006l ). For Haro 29, we used only U, B and J bands 
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from iHunter fc Elmegreenl (120061 ) because the V-band photometry is inconsistent with the 
others. For N GC 4861, we used ugiz data from SDSS. The giant clump in NGC 2366 was 
also studied by Kennicutt et al.l ( 1980 ). The giant clump in NGC 4861 is known as Mrk 59 
and was studied by llzotov et al.l (j2009bl ) and others. 



Although there is HI gas present in some of the clumps, and perhaps even molecular gas 
connected with current star formation, we do not include gas in the cl ump masses. In DDO 



155, for example, the HI m ass in the clump is a few times 10 5 M & (ICarignan et al.l 11990 



Begum fc Chengalurl 120031 ) . which is comparable to the stella r mass. In NGC 23 66, an HI 
cloud at the position of the clump contains several xlO 5 M (IHunter et al.ll200ll ). which is 
~ 20% of the stellar mass. Considering the possible addition of gas, the clump masses and 
mass fractions given in Table [T] are lower limits. 

As part of the fits for c lump mass, we also obtained crude star formation histories in the 
clumps (IZhang et al.l 120 111 ). These are determined as relative stellar masses younger than 
0.1 Gyr, in the time interval between 0.1 Gyr and 1 Gyr, and older than 1 Gyr. For the five 
galaxies, the relative masses in the intervals (< 0.1 Gyr, 0.1 — 1 Gyr, > 1 Gyr) are, Mrk 
178: (0.42, 0.49, 0.09), DDO 155: (0.11, 0.23, 0.66), Haro 29: (0.17, 0.15, 0.68), NGC 2366: 
(0.28, 0.5, 0.22), and NGC 4861: (0.27, 0.50, 0.23). Evidently the SEDs indicate significant 
clump components older than 1 Gyr even after background disk subtraction. The dominant 
appearance of these clumps in the J-band (Fig. 2) suggests the same thing. If these old 
massive components are really present, then they would have to be gravitationally bound to 
the clump and the clump would have to be long-lived. We note that gravitational instabilities 
in a disk of gas and stars can collect both gas and a significant mass of background field 
stars int o a clump when the velocity dispersions and densities of the two components are 
similar (lElmegreenl 1201 ll ) . Background field stars also fall into t he clump and get trap ped 
because of the changing gravitational potential as its mass grows (IFellhauer et al.l 120061 ) . In 
our sample, the clumps that are relatively closest to the center (in Mrk 178, DDO 155 and 
Haro 29) contain the highest fraction of old stars. This suggests a larger total age for the 
more centralized clumps than for the more peripheral clumps, which is consistent with a 
history of inward migration. 

The 9th column in Table [T] gives the ratio between the clump mass and the galaxy 
dynamical mass inside the clump radius. As shown in the previous sections, if this ratio is 
larger than a few percent, the clump could significantly perturb the surrounding disk and 
cold halo particles, leading to the loss of clump orbital angular momentum in less than ~ 1 
Gyr. The tabulated mass fractions are in this range. The timescale for their migration is in 
the last column, assuming NFW and Burkert profiles in two cases, and using the observed 
rotation speed at the clump position. The mass fractions are higher and the timescales are 
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smaller when the clumps are relatively close to the center, because only the inner parts of 
the galaxies are included in M dyn . 



3.2. Galaxy Thickness and Scale Length Ratios 

Clump accretion can thicken the central regions because the stellar orbital ene rgy gets 



mixed into three dimensions during the final merger phase (IBournaud et al.l 120071 ) . What 
is important is the ratio of the disk scale height, H = a 2 /(nGU), to the disk scale length 
i?d- Here, a is the perpendicular velocity dispersion in the central region of the BCD, and 
S is the central mass column density of the disk. For re ference values a = 10 km s~ l and 



£ = 10 M pc~ 2 , we obtain H = 740 pc. Most BCDs in Hunter & Elmegreenl (120061 1 have 
i?d ~ 500 pc or less, so H and R<± are comparable. This means the inner parts of BCDs are 
3D objects like a bulge. 

Detailed consideration of the BCDs in Table 1 confirm that the inner disk thicknesses 
are comparable to or larger than the inner disk scale lengths. Putting dimensions into the 
thickness equation, we get 

H = 740(^/10 km s _1 ) 2 (E/10 M pc^V- (9) 



Zhang et al.l ( 120111 ) determined disk stellar mass densities and scale lengths from SED fits. 
Values are given in Tabled The average H/R^ ~ 2.6, so the BCDs in this study should have 
relatively thick inner regions. The ratio would be larger for larger perpendicular velocity 
dispersions - the assumed value of 10 km s _1 for stars seems to be a lower limit. We note 
that the BCDs with giant clumps closest to the center have the highest ratios of height to 
length. 

Height-to-length ratios greater than unity in Table [2] are difficult to understand. They 
would be smaller if additional mass were in the disk. This suggests that some of the BCDs 
in our survey have a considerable mass column density of gas in the inner disk, perhaps 
comparable to or larger than the stellar column density. Such high masses of gas might 
be expected for the clumps in which the starbursts are occurring (e.g., larger than several 
hundred M pc -2 in molecules, which is typical for local giant molecular clouds), but there 
might also be a dense molecular and atomic intercloud medium where the average exceeds 
the average stellar value of 10 M pc~ 2 . Alternatively, a high filling factor of star-forming 
gas clumps that individually have mass column densities in excess of ~ 100 M pc -2 could 
produce an average gas column density in the inner part that exceeds the stellar column 
density. This could explain why the BCDs in Table[2]that have their massive clumps closest to 
the center also have the largest height-to- length ratios, i.e., these galaxies have higher S than 
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we assume because of contributions from molecular and dense atomic material in clumps. 
Massive clump accretion like that discussed here would drive significant gas accretion, not 
only in the clumps but also of the interstellar material between the clumps, which gets 
dragged along with the clumps by gravitational and magnetic forces. 

The mass column densities of inner disk HI gas have been obser yed for most of these 
galaxies. For both DDO 155 and NGC 2366, it is ~ 10 M Q pc" 2 flCarignan et~aD Il990 



Hunter et al.l 120011 ). Haro 29 has a hole in the ce ntral HI but the cl ump is very close and it 
has an average column density of ~ 20 M Q pc -2 (IStil fc Israeli 120021 ) . N GC 4861 has a larg e 
HI concentration in the center with a column density of ~ 30 M & pc -2 ( IThuan et al.l 12004 ) . 
Some of these values are larger than the corresponding central stellar mass column density 
by a factor 2 or more, which lowers H/R^ in proportion. Further studies of the gas column 
densities in the centers of BCD galaxies should clarify their relative thicknesses. 



4. Summary 

Young stellar clumps that form by gravitational instabilities in a galaxy disk can have 
such a high mass relative to the enclosed galaxy mass that they produce dynamically sig- 
nificant torques on the halo stars and cold dark matter particles, on the disk, and on each 
other. If the clump mass fraction exceeds a few percent, then these torques can drive an 
inflow of the clump's amount of mass in less than 1 Gyr. This process has been suggested 
for the formation of bulges in disk galaxies at high redshift, but it may apply also to local 
clumpy galaxies. Because of the general tendency for downsizing, in which active star for- 
mation occurs in galaxies with ever smaller masses as the universe ages, the clumpy phase 
now is mostly limited to dwarfs. We suggest that BCDs are an example of a local clumpy 
star-bursting galaxy in which the clumps are large enough to drive significant accretion in 
a Gyr or less. This would explain the dense stellar inner disks of these galaxies, and the 
prolonged star formation near the center. 

The BCDs in our sample also have relatively thick inner regions, reminiscent of bulges 
in spiral galaxies. They are even a little too thick if only the stellar surface densities are 
considered. This suggests there could be a dense atomic or molecular component in the inner 
region that has an average surface density comparable to or exceeding that from stars. 

This work was funded in part by the National Science Foundation through grants AST- 
0707563 and AST-0707426 to DAH and BGE. HZ was partly supported by NSF of China 
through grants #10425313, #10833006 and #10621303 to Professor Yu Gao. We are grateful 
to the referee for helpful comments. 
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Table 1. Sample BCD Galaxies and their Clump Properties a 



Galaxy 


D 


logM s 


logM b 


logM c 


r c 


Aperture 


log M dyn (r c ) 


M c /M dyn (r c ) 


In A 


T 




Mpc 


M Q 


M Q 


M e 


kpc 


kpc 


M© 






Gyr 


Mrk 178 


3.9 


7.04 


7.39 


5.13 


0.39 


0.32 


6.60 


0.035 


3.4-4.3 


0.50-1.4 


DDO 155 


2.2 


6.47 


7.22 


5.46 


0.21 


0.24 


6.78 


0.048 


3.2-3.9 


0.12-0.37 


Haro 29 


5.9 


7.16 


8.06 


6.33 


0.27 


0.74 


6.26 


1.17 


1.3-2.3 


0.03-0.06 


NGC 2366 


3.4 


7.84 


9.04 


6.23 


1.31 


0.93 


8.08 


0.014 


3.6-4.0 


1.3-4.2 


NGC 4861 


7.6 


8.04 


8.83 


6.89 


2.07 


0.67 


8.48 


0.026 


1.8-1.8 


1.7-6.3 



a D is the distance, M B is the galaxy stellar mass, M\, is the galaxy baryonic mass, M c is the clump stellar mass, r c is the 
clump galacticentric radius, Aperture is the aperture size used for clump photometry, M dyn is the galaxy dynamical mass inside 
r c , A is the Coulomb factor, and T is the clump accretion time. For the latter two, we assume £ = 0.2 and a rotation curve 
slope f) = 0.5 in the first case (NFW core), and £ = 0.03, /9 = 1 in the second case (Burkert core), with factors of 0.40 and 0.33 
in T, respectively, to account for the decrease in M dyn with radius. 
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Table 2. Inner Scale Heights and Lengths a 



Galaxy 


£ 


H 




H/R d 




M pc- 2 


kpc 


kpc 




Mrk 178 


4.9 


1.8 


0.27 


6.7 


DDO 155 


7.4 


1.0 


0.22 


4.5 


Haro 29 


39 


0.19 


0.20 


0.95 


NGC 2366 


0.66 


1.3 


3.7 


0.35 


NGC 4861 


10 


0.74 


1.0 


0.74 



is the inner disk scale height assuming a perpendicular velocity dispersion of 10 km s 1 and the observed stellar mass 
column density, S; R d is the inner disk scalelength. 
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Radial Position, x=r/r s 



Fig. 1. — Solutions to various parameters connected with the Burkert (1995) dark matter 
density profile, which has a constant-density core; £ is the dynamical friction parameter in 
the parentheses of equation [1], v is the rotation speed, a is the 3D velocity dispersion, /3 is 
the slope of the rotation curve, and T is the prefactor in equation |2j The dashed line shows 
(5 for a Navarro et al. (1996) profile. 



-27- 



Haro29 U+ 

_ • 

- (S) 

500 pc 


Mrkl78 U+ 

500 pc 


DD0155 U+ J 

* 

• * 

500 pc 


NGC2366 U+.I 

4 

500 pc 


| NGC4861 u+ 
500 pc '^H''' 


■ 



Fig. 2. — Five clumpy irregular galaxies used to study relative clump mass and possible 
central migration from tidal torques. The images are a combination of U and J-band from 
Hunter & Elmegreen (2006) except for NGC 4861, which is a combination of u and z-band 
from SDSS. The measured clumps are indicated by elliptical contours and the centers of the 
outer V-band isophotes are indicated by "x" . 



